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KNOTTED STRUCTURES IN HIGH-ENERGY BELTRAMI 
FIELDS ON THE TORUS AND THE SPHERE 

ALBERTO ENCISO, DANIEL PERALTA-SALAS, AND FRANCISCO TORRES DE LIZAUR 


Abstract. Let 5 be a finite union of (pairwise disjoint but possibly knotted 
and linked) closed curves and tubes in the round sphere S 3 or in the flat 
torus T 3 . In the case of the torus, S is further assumed to be contained in a 
contractible subset of T 3 . In this paper we show that for any sufficiently large 
odd integer A there exists a Beltrami field on § 3 or T 3 satisfying curl u = A u 
and with a collection of vortex lines and vortex tubes given by S, up to an 
ambient diffeomorphism. 


1. Introduction 

An incompressible fluid flow in R 3 is described by its velocity field u(x, t), which 
is a time-dependent vector field satisfying the Euler equations 

3 ( u+(m-V)u = -VP, divu = 0 

for some pressure function P(x,t). When the velocity field does not depend on 
time, the fluid is said to be stationary. This paper concerns stationary solutions of 
the Euler equations, which describe equilibrium configurations of the fluid. 

A central topic in topological fluid mechanics, which can be traced back to Lord 
Kelvin in the XIX century m, concerns the existence of knotted stream and vortex 
structures in stationary fluid flows. The most relevant of these structures are the 
stream lines, vortex lines and vortex tubes of the fluid. We recall that a stream line 
and a vortex line are simply a trajectory (or integral curve) of the velocity field u 
and the vorticity w := curlu, respectively, while a vortex tube is the interior do¬ 
main bounded by an invariant torus of the vorticity. The existence of topologically 
complicated stream and vortex lines is a central topic in the Lagrangian theory of 
turbulence and in magnetohydrodynamics, and has been studied extensively in the 
last decades (see e.g. mm for recent accounts of the subject). 

Our understanding of the set of stationary states of the Euler equations in three 
dimensions is much more limited than in the two-dimensional situation mm - in 
particular, the existence of stationary solutions in M 3 having stream lines, vor¬ 
tex lines and vortex tubes that are knotted and linked in arbitrarily complicated 
ways has been established only very recently mum- Following a suggestion of 
Arnold mm related to his celebrated structure theorem, to prove these results one 
does not consider just any kind of solutions to the stationary Euler equations but 
a very particular class that are called Beltrami fields. A Beltrami field in R 3 is a 
vector field satisfying the equation 

curl it = Xu 


( 1 . 1 ) 
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for some nonzero constant A. Notice that stream lines and vortex lines coincide 
in the case of a Beltrami field, and that a Beltrami field is automatically smooth 
(even real analytic) by the elliptic regularity theory. 

The stationary solutions in R 3 that one can construct using the techniques 
in mm fall off at infinity as l/|x|, this decay being sharp for Beltrami fields 
but not fast enough for the velocity to be in the energy space L 2 (R 3 ). In fact, 
the incompressibility condition ensures that there are no Beltrami fields in R 3 with 
finite energy even if the proportionality factor A is allowed to be nonconstant, as 
has been recently shown in [US- 

On the contrary, Beltrami fields in a closed Riemannian 3-manifold M (or a 
bounded domain of R 3 ) are stationary solutions to the Euler equations that do 
have finite energy. If S is a union of (possibly knotted and linked) closed curves 
and embedded tori in the 3-sphere, in this setting one can use contact topology to 
show m that there is a Riemannian metric g on the sphere with an associated 
Beltrami field u having a collection of vortex lines and vortex tubes given precisely 
by S. The main ideas of the proof are that the Reeb field of a contact form is in fact 
a Beltrami field in some adapted metric and that one can indeed construct contact 
forms on the sphere whose Reeb fields have the collection of periodic trajectories and 
invariant tori given by S. Notice that, as it is a Reeb vector field, a Beltrami field 
obtained in this fashion does not vanish. Conversely, any nonvanishing Beltrami 
field on the sphere is the Reeb vector field of some contact form, so in particular it 
must possess a closed vortex line [12]. 

Our goal in this paper is to establish the existence of knotted and linked vortex 
structures in Beltrami fields on compact manifolds with a fixed Riemannian metric. 
Specifically, we will consider Beltrami fields in the flat 3-torus T 3 and in the unit 3- 
sphere S 3 ; in fact, the former is the most fundamental space considered in the fluid 
mechanics literature other than R 3 and the latter is perhaps the simplest example 
of a closed Riemannian 3-manifold from a geometric point of view. 

It is worth emphasizing that, for a fixed Riemannian structure, the problem is 
much more rigid than when one can freely choose a metric adapted to the geom¬ 
etry of the set of lines and tubes that one aims to recover from the trajectories 
of a Beltrami field. An obvious reason is that, analytically, Beltrami fields in a 
closed Riemannian manifold arise as eigenfields of the curl operator, which defines 
a self-adjoint operator with discrete spectrum and a dense domain in the space 
of divergence-free L 2 fields. In the context of spectral theory, the proportionality 
constant A, or rather its absolute value, can be thought of as the energy of the 
Beltrami field, although of course it is in no way related to the L 2 norm of the 
latter. 

Our main theorem asserts that there are “many” Beltrami fields u in the sphere 
and in the torus with vortex lines and vortex tubes of any link type. Furthermore, 
these structures are structurally stable in the sense that any vector field on the torus 
or the sphere which is sufficiently close to u in the C 4 norm and which preserves 
some smooth volume measure will also have this collection of periodic trajectories 
and invariant tori, up to a diffeomorphism. To state this result precisely, let us 
call a tube the closure of a domain (in 8 3 or T 3 ) whose boundary is an embedded 
torus. Throughout, diffeomorphisms are of class C°°, curves are all assumed to be 
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non-self-intersecting, and we will agree to say that an integer is large when it is 
large in absolute value. 

Theorem 1.1. Let S be a finite union of (pairwise disjoint, but possibly knotted 
and linked) closed curves and tubes in § 3 or T 3 . In the case of the torus, we also 
assume thatS is contained in a contractible subset of T 3 . Then for any large enough 
odd integer A there exists a Beltrami field u satisfying the equation curlu = Xu and 
a diffeomorphism $ of § 3 or T 3 such that $(5) is a union of vortex lines and vortex 
tubes of u. Furthermore, this set is structurally stable. 

An important observation is that the proof of this theorem yields a reasonably 
complete understanding of the behavior of the diffeomorphism <h, which is, in par¬ 
ticular, connected with the identity. Oversimplifying a little, the effect of $ is to 
uniformly rescale a contractible subset of the manifold that contains S to have a 
diameter of order 1/|A|. In particular, the control that we have over the diffeo¬ 
morphism $ allows us to prove an analog of this result for quotients of the sphere 
by finite groups of isometries (lens spaces). Notice that $( 1 !?) is not guaranteed 
to contain all vortex lines and vortex tubes of the Beltrami field. It is also worth 
mentioning that, if S only consists of curves, the condition that the perturbation of 
the Beltrami field be volume-preserving is not necessary for the structural stability 
of <k(S), and the smallness in C 4 can be replaced by a C 1 condition. 

In § 3 and T 3 , Theorem 11.11 proves a conjecture of Arnold [Tj asserting that there 
should be Beltrami fields having stream lines with complicated topology. Further¬ 
more, it should be noticed that the helicity %(curl u) of the vorticity is proportional 
to its eigenvalue A so the Beltrami fields constructed in the main theorem have very 
large helicity. More precisely, the scale-invariant quantity %(curlu)/||u||| 2 , which 
is given by A in the case of a Beltrami field, becomes arbitrarily large. This is fully 
consistent with Moffatt’s interpretation [15] of helicity as a measure of the degree 
of knottedness of the vortex lines in the fluid flow. 

The proof of the theorem involves an interplay between rigid and flexible prop¬ 
erties of high-energy Beltrami fields. Indeed, rigidity appears because high-energy 
Beltrami fields in any 3-manifold behave, locally in sets of diameter 1/A, as Bel¬ 
trami fields in R 3 with parameter A = 1 do in balls of diameter 1. The catch here 
is that, in general, one cannot check whether a given Beltrami field in R 3 actually 
corresponds to a high-energy Beltrami field on the compact manifold. To prove 
a partial converse implication in this direction fTheorem 12. Ill , it is key to exploit 
some flexibility that arises in the problem as a consequence of the fact that large 
eigenvalues of the curl operator in the torus or in the sphere have increasingly high 
multiplicities. For this reason the proof does not work in a general Riemannian 
3-manifold. 

One should notice that the techniques introduced in [ 3 ® to prove the existence 
of Beltrami fields in R 3 with a prescribed set S of closed vortex lines and vortex 
tubes do not work for compact manifolds. The reason is that the proof is based 
on the construction of a local Beltrami field in a neighborhood of S, which is 
then approximated by a global Beltrami field in R 3 using a Runge-type global 
approximation theorem. For compact manifolds the complement of the set S is 
precompact, so we cannot apply the global approximation theorem obtained in [7] 
®. In fact, as is well known, this is not just a technical issue, but a fundamental 
obstruction in any approximation theorem of this sort. This invalidates the whole 


4 ALBERTO ENCISO, DANIEL PERALTA-SALAS, AND FRANCISCO TORRES DE LIZAUR 


strategy followed in 13 El and makes it apparent that new tools are needed to prove 
the existence of Beltrami fields with geometrically complex vortex lines and vortex 
tubes in compact manifolds. 

The paper is organized as follows. In Section [2] we will prove the main theorem 
assuming that Theorem 12. II holds. Theorem 12.II will be proved in Section [3] in the 
case of the sphere, with the proof of some technical results relegated to Sections0] [HI 
and in Section [7] in the case of the torus. The paper concludes with some remarks 
that we present in Section [51 where in particular we prove an analog of the main 
theorem for lens spaces. 


2. Proof of the main theorem 


For the ease of notation, we shall write M 3 to denote either T 3 (the standard 
flat 3-torus, (R/27rZ) 3 ) or § 3 (the unit sphere in R 4 ). A Beltrami field u in M 3 is 
an eigenfield of the curl operator, which satisfies 


curl u = Xu , 


for some nonzero constant A. It is well-known that the spectrum of curl in the 
sphere are the integers of absolute value greater than or equal to 2, while in T 3 
consists of the real numbers of the form 


A = ±|fc| 


for some k £ Z 3 . In particular, the spectrum of curl in T 3 contains the set of 
integers. Here and in what follows, | • | denotes the usual Euclidean norm of a 
vector. 

The following theorem, whose proof is presented in Section [3l shows that a Bel¬ 
trami field v in R 3 can be approximated, up to a suitable rescaling, by a high-energy 
Beltrami field u in M 3 . This fact is key to the proof of Theorem 11.11 as it implies 
that the dynamics of any Beltrami field of R 3 in compact sets can be reproduced 
in a small ball of M 3 by a high-energy Beltrami field on the manifold, provided 
that the dynamical properties under consideration are robust under suitably small 
perturbations. For concreteness, we will henceforth assume that A is positive; the 
case of negative A is completely analogous. 

For the precise statement of the theorem, let us fix an arbitrary point po G M 3 
and take a patch of normal geodesic coordinates T : B -» B centered at pq. Here 
and in what follows, B p (resp. B p ) denotes the ball in R 3 (resp. the geodesic ball 
in M 3 ) centered at the origin (resp. at po) and of radius p, and we shall drop the 
subscript when p = 1. The theorem will be then stated in terms of the vector 
field ’F*it on H, which is just the expression of the Beltrami field u in local normal 
coordinates. If u k (x ) are the three components of 'F*t( in the Cartesian basis {ej}? =1 


of R 3 , i.e. 


3 



4=1 

we will make use of the rescaled vector field 
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Theorem 2.1. Let v be a Beltrami field in R 3 , satisfying curli; = v. Let us fix 
any positive numbers e and m. Then for any large enough odd integer X there is a 
Beltrami field u, satisfying curlrt = Xu in M 3 , such that 


( 2 . 1 ) 



< e. 

C m (B) 


Let us now show how this result can be exploited to prove the main theorem. 
For this, let be a diffeomorphism of M 3 mapping the set S into the ball Bi/a> 
and the ball 1m into itself. (In S 3 , the existence of such a diffeomorphism is 
trivial, while in the case of T 3 it follows from the assumption that S is contained 
in a contractible set.) Furthermore, given a positive number A let us denote the 
rescaling with factor A by ©a(^) := Ax. We can now define a set S' of finitely 
many closed curves and tubes in the ball B as 

: = (0 A o^o $')(£). 

The following result is a straightforward consequence of the main theorem in [§] : 

Theorem 2.2. There is a Beltrami field v in R 3 satisfying curlp = v and an 
orientation-preserving diffeomorphism $0 0 /R 3 , which coincides with the identity 
in the complement of B, such that $0 {S') is a union of vortex lines and vortex tubes 
of v. Furthermore, this set is structurally stable. 


Proof. It was shown in [9] that there is a Beltrami field v in R 3 , satisfying 

curie = Xv 

for some small positive constant A < 1, and an orientation-preserving diffeomor¬ 
phism $ of R 3 that is the identity in the complement of B such that $(5') is a 
set of closed vortex lines and vortex tubes of v. The closed vortex lines are elliptic 
trajectories of v and the boundaries of the vortex tubes are KAM-nondegenerate 
invariant tori of v. The theorem follows setting v(x ) := v(x/X), which satisfies the 
equation curlu = v in R 3 , and noticing that (0 A o <F)(S') is a set of closed vortex 
lines and vortex tubes of v. Since this set is contained in B because A < 1, it is 
standard that there exists a diffeomorphism $0 of R 3 mapping S' onto 0^ o < h(5') 
which is the identity in the complement of B. The closed vortex lines in the set 
, I>o(lS , ) are structurally stable under C 1 -small perturbations by the elliptic per¬ 
manence theorem, while the vortex tubes are structurally stable under C' 4 -small 
volume-preserving perturbations by the KAM theorem. □ 


Let us now combine Theorems I2A1 and !2.2l to conclude the proof of Theorem ll.il 
Theorem 12.11 guarantees that, for any large enough odd integer A, the Beltrami 
field v constructed in Theorem 12721 can be approximated in the sense of Eq. (E3D 
by a Beltrami field u defined on M 3 . Then it is not hard to see that the structural 
stability of the set $o(«S0 of closed vortex lines and vortex tubes of v implies the 
existence of a diffeomorphism $1 of R 3 , which is the identity in the complement 
of B, such that C B is a set of structurally stable closed vortex lines and 

vortex tubes of the rescaled held 

*•“(1)• 


(2.2) 
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Indeed, because of the elliptic permanence theorem, this claim is immediate in the 
case of closed vortex lines provided that the number m appearing in the approxi¬ 
mation estimate m is at least 1. For the case of vortex tubes one can use that 
the Beltrami field u is divergence-free in M 3 , which ensures that the field (12.21) pre¬ 
serves a smooth volume 3-form in B that is a small perturbation of the Euclidean 
one, namely 

($*M) ("0 = Mo + 0(X 1 ). 

Here p, and p o respectively denote the canonical volume 3-forms of M 3 and R 3 . 
Hence, taking m ^ 4 in the approximation estimate m, this enables us to apply 
the KAM theorem for volume-preserving fields in R 3 , which ensures the existence of 
the aforementioned diffeomorphism yielding the desired set of vortex tubes of the 
rescaled field (12.21) . (For the benefit of the reader let us recall that, in order to prove 
this KAM result, one takes a Poincare section transversal to the tube of v under 
consideration, thereby reducing the problem to perturbations of a nondegenerate 
twist map of the annulus with the intersection property. It is then standard that 
one can apply a Moser-type twist theorem to guarantee the preservation of the 
invariant tori. The details, which go as in [9J Section 7.4], are omitted.) 

It follows from the above discussion that the diffeomorphism $ of M 3 can be 
then defined as 


$(x) 


$'(*) if x£ fc'-HBi/*), 

(4/ _1 o o o 0 A o o $')(a;) if x £ . 


The set <f>(<S) is then the union of structurally stable closed vortex lines and vortex 
tubes of the Beltrami held u, so the main theorem follows. 


3. Proof of Theorem 12.II in the sphere 

In this section we show that for any Beltrami held v in R 3 there exists a Beltrami 
held u in § 3 satisfying curl it = Xu whose dynamics in a ball of radius A -1 is very 
close to the dynamics of v in the unit ball. The proof is divided in three steps. In 
the first step we show that the Beltrami held v can be approximated in B by a held 
w that is a hnite sum of spherical Bessel functions jo(\x — x n \) centered at different 
points x n £ R 3 ('Proposition 13.11) . The held w is not generally a Beltrami held, 
however. In the second step we show that one can take three spherical harmonics 
Y lt Y 2 ,Y 3 in § 3 of energy A(A—2) whose behaviors inaball ofradius 1/A respectively 
correspond to those of the three components of the held w in a ball of radius 1, 
provided that A is large enough (Proposition 13.21) . Finally, in the third step we 
construct a Beltrami held u in § 3 of energy A, using as key ingredients the spherical 
harmonics Yfc and a basis of Hopf helds, so that u approximates the held v in the 
sense of Eq. (18.11) (Proposition 13.31) . 

For notational convenience, in this section we will write A := A — 2. Notice that 
A is then a large integer. 

Step 1: Approximating the Beltrami field v by sums of shifted spherical Bessel 
functions. The hrst step of the proof of Theorem 12 .1 1 consists in showing that there 
is a hnite sum w of spherical Bessel functions jo centered at different points that 
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approximates the Beltrami field v in the unit ball of R 3 . The field w is not a 
Beltrami field but, just as v, it satisfies the Helmholtz equation 

Aw + w = 0 . 


Proposition 3.1. For any 8 > 0, there is a finite radius R and finitely many 
constants {c n }^ =1 C R 3 and {x n }^ =1 C Bn such that the field 

N 

W := ^2 c njo{\x - x n \) 
n—1 

approximates the Beltrami field v in the ball B as 


\\v - w\\ c ™+2( B ) < 6 . 


The proof of this proposition will be presented in Section [4] 


Step 2: Approximating the field w by high-energy spherical harmonics. Let us write 
the vector field w in terms of its components w l in the Cartesian basis {e ;} 3 _ 1 of R 3 : 

3 

w = w 'ei. 

i—l 

Each component w l is a solution of the Helmholtz equation Aro' + w' = 0 in R 3 . We 
now show that for any large enough integer A, there exists a spherical harmonic Yj 
on § 3 with energy A(A + 2) that behaves in the ball Bi/a as w l does in the unit ball. 
The proof of this result is based on the asymptotic expressions for the ultraspherical 
polynomials, which are the building blocks for any spherical harmonic on § 3 , and 
exploits in a crucial way the expression for w as a finite sum of spherical Bessel 
functions that we obtained in Step 1: 


Proposition 3.2. Given any positive constant 8, for any large enough integer A 
there is a spherical harmonic Y* on § 3 with energy A(A + 2) such that 


w i -1 

Yi 0 T -1 f—) 


VA / 


C m + 2 (B) 


< 8 . 


The proof of this proposition is given in Section [5] 


Step 3: Construction of the Beltrami field on § 3 using spherical harmonics and 
Hopf fields. Let us consider the three positively oriented orthonormal Hopf vector 
fields in § 3 that, in terms of the Cartesian coordinates of R 4 , are explicitly given 

by 

hi := (—x 4 , 23 , -x 2l xi ), 
h 2 := (—£3, —x 4 , X\,x 2 ), 
h 3 := (-x 2 ,x 1 ,x 4 ,-x 3 ). 

It is well known that they are curl eigenfields with eigenvalue 2, that is, 

curl hj = 2 hi. 

We have taken the the Cartesian basis ei of R 3 so that ’L*hi(0) = e 

In the following proposition we show how to construct a Beltrami field on § 3 
using the spherical harmonics Yi obtained in Proposition 13.21 and the Hopf fields hi 
so that it approximates the Beltrami field v in a suitable sense. 
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Proposition 3.3. The vector field on the sphere 

u := —^ curl(curl +A) (Yihi + Y 2 h 2 + P 3 / 13 ) 
is a Beltrami field satisfying curlu = (A + 2 )u and approximates v as 

provided that A is sufficiently large. 

Here C is a constant depending on m but not on S. Since rescaling by A is 
essentially equivalent to rescaling it by A because 


< CS, 


C™(R\ 


11 
A ~ A 



Theorem [57T] then follows from Proposition l3. 31 provided A is sufficiently large and 6 
is chosen small enough for CS not to be larger than e/2. The proof of Proposition l3.3l 
is given in Section [G] 


4. Proof of Proposition 13.II 


Since the Beltrami field v satisfies the Helmholtz equation Av + v = 0, upon 
expanding the components of v in a series of spherical harmonics it is elementary 
to realize that v can be written in the ball B 2 as a Fourier-Bessel series of the form 

OO l 

(4.D -EE Sim jl (t) P/m (a^) 

1=0 m=—l 


that converges in L 2 (B 2 ). Here r := \x\ £ R + and w := x/r £ S 2 are spherical 
coordinates, ji is the spherical Bessel function, Yi m are the spherical harmonics and 
bi m £ R 3 are constant vectors. 

Since the series 63D converges in L 2 {B 2 ), for any S' there is an integer Iq such 
that the finite sum 

lo l 

Vi := E E him ji[r)Y lm {uj) 

1=0 m=—l 

approximates the field v in an L 2 sense, that is, 

(4.2) ||ffi-u||L2 (b 2 )<6'. 


Next, let us observe that the properties of the spherical Bessel functions imply 
that the field Vi falls off at infinity as |vi(a;)| < C/ |x|. In particular, it then follows 
from Herglotz’s theorem (see e.g. 13] Theorem 7.1.27]) that v\ can be written as 
the Fourier transform of a distribution supported on the unit sphere of the form 

(4.3) Vl (x)=[ h(f)e lx< dcr(f ), 

ds 2 

where da is the area measure induced on the unit sphere S 2 := {(el 3 : |£| = 1} 
and /i is an Revalued function in L 2 (§ 2 ). 

By the density of smooth functions in L 2 (§ 2 ), we can approximate fi by a smooth 
function f 2 : S 2 —> R 3 so that their difference is bounded as 

ll/i ~ /2 ||l 2 (§ 2 ) < S'. 
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Therefore the field 

(4.4) v 2 (x) := f f 2 (0 e lx< da (£), 

Js, 2 

approximates V\ uniformly, as for any x £ R 3 the Cauchy-Schwarz inequality yields 


(4.5) \v 2 (x) - V!(x)\ = 


Is 2 


(f2(0~fme ix< da(0 


< C\\f 2 - AIIl^s 2 ) < CS'. 


Our next objective is to show that for any S' there is a radius R > 0 and finitely 
many constants { c n }^ =1 C R 3 and {x n }^ =1 C Br such that the restriction to the 
unit sphere of the smooth field in R 3 


N 


m - V ■ 


0 —ix n ^ 


n —1 


approximates the field f 2 in the C° norm, that is, 

(4-6) 11/- /2||c°(S 2 ) < S'. 

To prove this claim, we first extend f 2 to a smooth vector field g : M 3 —> R 3 with 
compact support. This can be done setting 

g{0 ~ X(l^l)/2^||| 

with y(s) a smooth function which is equal to 1 , say, if |s — 11 < \ and vanishes 
for | s — 11 > 4. Since the Fourier transform ^ of <7 is Schwartz, we easily infer that 
there is a radius R such that the L 1 norm of <7 is essentially contained in Br in the 
sense that 

f \g(x)\dx<6'. 

It then follows that the Fourier integral representation of g can be essentially trun¬ 
cated to an integral over the ball Br , i.e., one has the uniform bound 


(4.7) 


sup 

£<EK 3 


5(0 - / 9(x) 

Jb r 




dx 


< 5’. 


Now, an easy continuity argument allows us to uniformly approximate the inte¬ 
gral 


[ g(x ) e~ ix < dx 
J Br 


by a finite sum 
(4.8) 


N 


m ■■= E 


Cn C 


— iXn‘£ 


with constants c n £ R 3 and x n £ Br in such a way that the error introduced in 
the approximation is bounded by 


(4.9) 


sup 


f g(x)e lx< dx-f(£) 
Jb r 


< S'. 


Indeed, let us cover the ball Br by finitely many closed sets {[/„}^ =1 with 
piecewise smooth boundaries and pairwise disjoint interiors such that the diameter 
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of each set is at most 5". The function e “'£<?( x) being smooth, it then follows 
that for each x,y € U n one has 

sup \g(x) e~ lx< - g{y) e~ ly< \ < CS" , 

{es 2 

where the constant C depends on <7 (and therefore on S') but not on 6". It is then 
straightforward that if x n is any point in U n and we set c n := g(x n ) \U n \ in (14.81) . 
one has 


sup 

feS 2 



N r 

<V/ sup \g(x) e~ ix '^ 


^ CS" 


g(x n )e dx 


where again C depends on S' and R but not on 5" or N. Hence one can take S " 
small enough so that CS" < S', thereby proving the estimate (1T91) . 

Putting together the estimates (1X71) and (14.91) we infer that 

ll/-<?llc°(s 2 ) <C<5', 

with a constant independent of S'. Since the restriction to § 2 of the function g is 
precisely fi, the estimate (14.61) then follows. 

Finally, if we define the vector field 

w(x) := [ /(£) 

we conclude from Eq. (14.61) that 

l|w-i>2||c°(H3) < [ 1/(0 - / 2 (C)! <M£) < CS', 

J s 2 

so we readily infer from Eqs. (14.21) and (14.51) the L 2 bound 
(4 - 10) 

||u — w 1 1 l 2 ^ C\\w — r>21|c°(R 3 ) + CJ\\v2 — Vi ||c°(m 3 ) + Iki — 0 Il 2 (b 2 ) < CO ■ 

Furthermore, as the Fourier transform of w is supported on S , w satisfies the 
Helmholtz equation 

Aw + w = 0 

in K 3 . Since the Beltrami field v also satisfies the Helmholtz equation Av + v = 0, 
standard elliptic estimates enable us to promote the L 2 bound (14.101) to the C m+2 
estimate 

\\v - m’IIc-+ 2 (b) < C\\v - w\\ L 2 (B2) < CS’, 
so the proposition follows upon choosing C5' < 5. 


N 

e ix < daft) = J2 C ' 

n= 1 


r . N 

/ e *(* X C< date) = y2c n jo(\x - x n \) 

J ® 2 n= 1 


5. Proof of Proposition 13.21 


For any positive integer A, let C A {t) be the ultraspherical (also called Gegen- 
bauer) polynomial of dimension 4 and degree A, which can be defined in terms of 
the Jacobi polynomials ^ as 


(5.1) 


C A (t) := 


yX r(A + l) ( 1 , 
2 F(A + §) A 


(*), 


where we are using the normalization (7a(1) = 1 for all A. 
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If p, q £ § 3 are two points in the 3-sphere, understood as the subset {|p| = 1} 
of R 4 , the addition theorem for ultraspherical polynomials shows that C'a (p ■ q) can 
be written as a linear combination of spherical harmonics. Specifically, 

2 (A+i) 2 

(5-2) C A (p ■ q ) = " + 1)2 Ya o(p) Y ^M) > 

where {Iaj}j =^" 1 ' 1 is an arbitrary orthonormal basis of spherical harmonics of en¬ 
ergy A(A + 2) and p ■ q denotes the scalar product in R 4 of the unit vectors p , q. 
Notice that (A 1 ) 2 is precisely the multiplicity of the eigenvalue A(A + 2) of the 
Laplacian on § 3 (or, equivalently, the dimension of the space of spherical harmonics 
of energy A(A + 2)). 

Let us write the i th Cartesian component of the vector field w as 

N 

W \ x ) = ^2 C nJo{\ x - X n \) , 

71=1 


where c l n is the i th component of the constant c n £ 
contained in the ball Br. Let us set, for any p £ § 3 , 

N 

Yi(p) ■= ^2 c « Ca ’ 


and the points x n are 


with 




Note that p n is well defined provided A is bigger than R. It is obvious from Eq. (15.21) 
that Yj is a spherical harmonic of energy A(A + 2). 

In order to study the asymptotic properties of the spherical harmonic Yi we first 
observe that, if we restrict our attention to points of the form 


p := 

with x £ Br and A > R 1 we then have 
(5.3) P- Pn = cos (distga (p,p n )) = cos 


\x - x n \ + 0(A 4 ) 
A 


as A —> oo. Here dists 3 (p,p„) denotes the distance between the points p and p n as 
measured on the sphere § 3 and the last equality stems from the fact that : IB —A? 
is a patch of normal geodesic coordinates. We will henceforth use the notation 


(5.4) 


Yi(x) := Y o $r 


Since for A large we have the asymptotic behavior 


r(A + i) _ 
r (A +1) " 

we conclude from Eq. (15.31) that 
(5.5) C A (p ■ Pn) = + 0(A“i) 


+ 0(A ~ 


»), 



x n \ + 0 (A 4 ) ^ 
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Now Darboux’s asymptotic formula for the Jacobi polynomials [T5J Theorem 8.1.1] 
implies 

-j=P ( k 2 \co S -)=^=j 0 (t) + °(A 1 ), 

which holds uniformly for compact sets (e.g., for |t| ^ 2 R). Therefore Eq. (15.41) can 
be written by virtue of Eq. (1531) as 


N 


n— 1 
N 



r \x - x n \ + 0 (A 3 )^ 


A )) 


= ^2c l n j 0 (\x-x n \)+O(A 1 ), 


n—1 


provided that A is sufficiently large and x,x n € Br. This proves that for any S' > 0 
and all A large enough we have the uniform bound 

(5.6) lk-^||c°(B) <#. 


To get the C m+2 bound stated in the proposition, we notice that the eigenvalue 
equation 

A Yi + A(A + 2 )Yi = 0 

for the spherical harmonic Y., in S 3 can be written in terms of the rescaled function Y) 
as 

A 0 h) + Yi = —AY), 

where the coordinates x are assumed to take values in B, A 0 := £3 d 2 . is the flat 
space Laplacian acting on the x coordinates and A is a scalar second-order operator 
of the form 

AY, := -2Yi + <?! DYi + G 2 D 2 ?i . 

Here the functions Gi(x, A) are (possibly matrix-valued) functions that depend 
smoothly on all their variables and whose derivatives are bounded independently 
of A for x € B, i.e., 

(5.7) sup | D^Gi(x, A)| < C a . 

XdzB 

Here the constant C a depends on the multiindex a but not on A. 

By construction, the function w l satisfies the Helmholtz equation 

A 0 w l + w l = 0 , 

and hence the difference w l — Yi satisfies the equation 

A 0 (w* - Y-) + (w i - Yi) = i A% . 

Therefore, in view of the uniform bounds (15.61) and (15.71) , standard elliptic estimates 
yield 

Ik - YiW C^(B) < C,\\w' - Yi\\ C0{B) + ^||Ayi|| C m,a (B) 
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which implies that 

II i T/- II / /'"'fsp I ^|| IPl II C' m T 2, ° ! . C 

||w ~ ^i||c m + 2 (B) ^ uo H- — - < 0 

provided that A is large enough (which in turn implies that S' is small). This 
completes the proof of the proposition. 


6 . Proof of Proposition 13.31 

We start by defining a vector field u on § 3 using the Hopf fields hi as 

u := Yi hi + Y 2 /12 + Y 3 ft. 3 , 

where the functions 1) are the spherical harmonics obtained in Proposition 13.21 In 
what follows it is convenient to work with differential forms, so let us denote by 
/3 and ai the 1-forms that are dual to u and hi, respectively, with respect to the 
canonical metric on § 3 . We recall that the dual of curl u is the 1 -form *d/3, with 
* being the Hodge star operator. 

In the following lemma we compute the action of the Hodge Laplacian on the 
1-form /3 using the properties of the Hopf fields: 

Lemma 6.1. The Hodge Laplacian of the 1-form /3 dual to u is 

—A/3 = A(A + 2) /3 + 2 * d/3, 

Proof. The 1-form /3 is given by /3 = Yj a.j , where summation over repeated indices 
is understood throughout. The Laplacian of /3 is then 

—A/3 := dd*l 3 + d*d/3 = —d * d * (1) on) + *d * d(Yj a,). 

Using that *doti = ‘Ion because on is the dual 1-form of the Hopf field hi, and 
that the differential of 1 ) can be written as dl) = hj(Yi) ctj, where hj(Yk ) denotes 
the action of the vector field hj on the scalar function Yj., we readily obtain 

d * d* (Yi a,;) = —d * (hj (Yi) aj A doti). 

Observing that ay A da.i = 2ay A *aj = 25jk /i, where p stands for the Riemannian 
volume 3-form on S 3 , it follows that 

(6.1) d* d* (YjOti) = d(hi(Yi)) = hjh^Y^aj . 

Analogously, a straightforward computation using that *(a 3 -Aai) = £juan, where 
£ju stands for the Levi-Civita permutation symbol, and the identity £i m i£jki = 
Sj'jS r nk SikSmj yields 

(6.2) * d(Yi af) = £j l ihj(Y i ) ai + 2Yj , 

(6.3) a d a d(Yi oti) — hj hj (Y)) oti T hihj (Y)) ay T 4 £ju hj (Yi) on - f- 4Y) oci. 
Finally, adding Eqs. m and m we obtain 

A/3 hjhi(Yi) cXj T hihj(Yi) ay hjhj(Yi) a% T 4c juhj (Yj) ai T 4Y) a^ 

= A(A + 2)Y) ai + 2£jnhj(Yi) ai + 4Y^ a^, 
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where we have used that AY) = — A(A + 2)Yi and that the commutator of Hopf 
fields is [hi,hj] = —2sijihi. The lemma then follows upon noticing that 

2Ejuhj(Yi) ai + 4 Yi on = 2 * df3 

by Eq. (16.21) . □ 

Using this lemma, it is easy to check that 

u := curl(curl+A)zt 

is a Beltrami field with eigenvalue A + 2. Indeed, if /3 is the dual 1-form of it, we 
obtain 

*d/3 = ^T A)/3 = —A -f A* d)/3 

= 2^2 * + A)/3 = (A + 2)/?. 

To prove the C m estimate of the proposition, it is convenient to introduce an 
auxiliary vector field in the unit ball B of R 3 as 

u(x) := Yi(x) e\ + Y 2 (x) e 2 + Y 3 (x) e 3 , 

where x £ B and Yi was defined in (15.41) . There is no loss of generality in choosing 
the orthonormal basis e; of M 3 compatible with the Hopf fields hi in the sense that 
4'*(/ii)(0) = ei . It is then easy to check that for x £ B one has: 

= B +x“' 

4'*(curlu)^—^ = A^curio u + + ^-Duj , 

/ • \ / G G G \ 

T*(curl curl J = A 2 ^ curio curio u + + -^-D 2 uj . 


Here curio denotes the Euclidean curl operator, acting on the variables x, and the 
functions Gi(x, A) are (possibly matrix-valued) functions that depend smoothly on 
all their variables and whose derivatives are uniformly bounded as 

(6.4) sup \D%Gi(x, A)| < C a . 

x(zB 


Here the constant C a depends on the multiindex a but not on A. 

These identities and the fact that (curio curio + curio)?; = 2v then permits us to 
write 


4 '*“(a) 


— V 


€ 


C m (B) 


-(curio curio + curio)(it — v) 


C, 


< C\\u - v\\ C-+ 2 (S) + -7-||« - w\\c™+*(B) 


C m (B) + V lU]lCm+2 ^ 

c „ 

A 1 

c c 

(6.5) + —Hu - w\\ C m+ 2 {B) + — \\v\\ C m+ 2 (B) . 

To conclude, notice that it stems from Propositions 13. II and 13.21 that 

II?; - it;||c m + 2 (B) < <5 
||u- w;||c”*+2 ( b) < 3<5, 
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so in particular 

||« - v\\c™+ 2 (B) < ||« - w\\c™+ 2 (B) + \\v - HI C™+ 2 (B) < 4(5 . 

Hence the proposition follows from the estimate (16.51) upon noticing that v is a 
fixed vector field (so its norm is independent of A) and choosing A large enough, 
which also allows us to take S as small as one wishes. 

7. Proof of Theorem 12.fl in the torus 

Arguing as in the proof of Proposition ^. II we can readily show that for any S > 0, 
there exists a vector field v± on M 3 that approximates the Beltrami field v in the 
ball B as 

(7.1) IK - Hc°(b) < H 

and that can be represented as the Fourier transform of a distribution supported 
on the unit sphere of the form 

vi(x)= [ f(0e* x da(0. 

J s 2 

Again § 2 denotes the unit sphere {£ € R 3 : |£| = 1} and / is a smooth Revalued 
function on § 2 . 

Let us now cover the sphere § 2 by finitely many closed sets {U n }^ =1 with piece- 
wise smooth boundaries and pairwise disjoint interiors such that the diameter of 
each set is at most 5'. We can then repeat the argument used in the proof of 
Proposition 13.II to infer that, if is any point in U n and we set 

Cn ■= f(£n) \U„\ , 

the field 

N 

w(x) := y^c n e <n ' x 

n= 1 

approximates the field vi uniformly with an error proportional to S': 

IH - II C°(b) < C5'. 

The constant C depends on S but not on S', so one can choose the maximal diam¬ 
eter S' small enough so that 

(7.2) |k-Ui|| C o(B) < (5. 

In turn, the uniform estimate 

Ik - Hc°(b) < Ik - ihllc°(B) + k - II G°(b) < 2 S 
can be readily promoted to the C m+ 2 bound 

(7.3) \\w-v\\ C rn.+ 2 {B ) <C5 . 

This follows from standard elliptic estimates as both w (whose Fourier transform 
is supported on S 2 ) and v satisfy the Helmholtz equation: 

Av + v = 0 , Aw + w = 0. 

Furthermore, replacing w by its real part if necessary, we can safely assume that 
the field w is real-valued. 
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Let us now observe that for any large enough odd integer A one can choose the 
points £n £ U n C S 2 so that they have rational components (i.e., £„ £ Q 3 ) and the 
rescalings A£ n are actually integer vectors (i.e., A£ n £ Z 3 ). This is because rational 
points £ £ § 2 fl Q 3 with A£ £ Z 3 are uniformly distributed on the unit sphere as 
A —> oo through odd values [ 6 ]. 

Choosing £ n as above, we are now ready to prove Theorem 12.11 in the torus. 
Without loss of generality, we will take the origin as the base point p, so that we 
can identify the ball B with B through the canonical 27r-periodic coordinates on the 
torus. In particular, the diffeomorphism ik : B — > B that appears in the statement 
of Theorem I2T1 can be understood to be the identity. 

Since A£„ £ Z 3 , it follows that the vector field 

N 

u{x) := ^2 c n e lAU ' x 

n =1 

is 27r-periodic (that is, invariant under the translation x — > x + 27r a for any vector 
a £ Z 3 ). Therefore it descends to a well-defined vector field on the torus T 3 := 
R 3 /(27rZ) 3 , which we will still denote by u. 

Since the Fourier transform of u if now supported on the sphere of radius A, u 
then satisfies the Helmholtz equation on the flat torus T 3 with energy A 2 , 

A u + A 2 u = 0 . 


A straightforward calculation then reveals that the vector field on the torus 

curl curl u + A curl u 


u := 


2A 2 


satisfies the equation 

curl u = A u , 

so it is a Beltrami field on T 3 with eigenvalue A := A. 
Let us now notice that, with some abuse of notation, 


A 


= w(x) 


for all points x, say, in the ball B. In particular, as the derivatives of the rescaled 
vector field u(-/ A) behave as 

curlrt^—^ =A curl w, 
curl curl u ^ = A 2 curl curl w , 

it then follows that 


( • 

\ 


A 2 curl curl w + A 2 curl w 

m (a 


C m (B) 

2A 2 t 


curlcurl(u; — v) + curl(u> — v) 


C m (B) 


C rn (B) 
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where we have used the identity curl curl v + curl v = 2v to pass to the second 
equality and the estimate (17.Ill) to derive the last inequality. The theorem then 
follows provided that <5 is chosen small enough for CS < e. 

8. Concluding remarks 

To conclude, let us make a few simple observations about our main result that 
follow from its proof: 

There are many Beltrami fields with closed vortex lines and tubes of a given link 
type. Indeed, since our construction works for any large enough odd integer A and 
Beltrami fields corresponding to different eigenvalues are L 2 orthogonal, there are 
many non-proportional Beltrami fields with closed vortex lines and tubes realizing 
any given link. 

In the sphere, the result holds true for any large enough eigenvalue A. Indeed, the 
fact that A is odd was never used in the proof of Theorem 12.II in § 3 (cf. Section [3]), 
so it stems that, given any finite union of closed curves and tubes S, for any 
integer A with |A| greater than certain constant Ao(«S) there is a Beltrami field 
with eigenvalue A having a structurally stable set of vortex lines and vortex tubes 
diffeomorphic to S. 

In our Beltrami fields on the sphere, knots and links appear in pairs. In fact, using 
the Hopf basis introduced in Section [3] any Beltrami field u on § 3 with 

eigenvalue A := A + 2, with A a nonnegative integer, can be written as 

U = Fi hi + i<2 /l2 + F3 /13 , 

where F t are smooth functions on the sphere. It is then easy to check using Eq. m 
that Fi must be a spherical harmonic of energy A(A + 2). Since such a spherical 
harmonic is known to have parity (—1) A , in the sense that 

F t (-p) = (~1) A F t (p) 

for all points p in the unit sphere § 3 , and the Hopf fields hi are odd (i.e., hi(—p) = 
~hi(p)), we conclude that a Beltrami field on the sphere with eigenvalue A has 
parity (—1) A+1 , so it is either even or odd. Therefore, the fact that $(5) is a 
set of vortex lines and vortex tubes of the Beltrami field u diffeomorphic to S and 
contained in a ball of small radius 1/A automatically implies that so is the antipodal 
set —$(5). 

The result carries over to lens spaces. In order to see why, the key is that in the 
sphere the statement of Theorem 12.11 can be refined to include localizations around 
different points of the sphere. More precisely, let us fix l points Pi,..Pi in § 3 , 
none of which are antipodal to another (that is, Pj / —Pk), and denote by T, : 
M(Pj,R 0 ) —>■ Br 0 a patch of normal geodesic coordinates centered at the point Pj. 
Here B (Pj,R 0 ) denotes the geodesic ball in the sphere of center Pj and radius 

R o '■= 7 ; mindist§3 (. Pj,P k ) • 

2 j^k 

The approximation theorem can then be stated as follows: 
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Theorem 8.1. Let {v>j}j—i be Beltrami fields in R 3 , satisfying curing = Vj . Let 
us fix any positive numbers e and m. Then for any large enough integer A there is 
a Beltrami field u, satisfying curl it = Ait in § , such that 


( 8 . 1 ) 

for all 1 ^ j ^ l. 



< e 

C m {B) 


Proof. Arguing as in Proposition 13.21 we infer that for any large enough integer A 
there are spherical harmonics Yjj of energy A(A + 2) such that 

*-'(-) 

where Wj is a vector field of the form 

N 

Wj = 'y ] Cj n jo(\ x ~ x jn |) 

71—1 


<S, 


C m + 2 (R) 


that approximates the Beltrami field Vj in C m+2 (B) as in Proposition 13.11 and 
Wj (1 ^ i ^ 3) denotes its i th Cartesian component. Noticing that the Jacobi 
polynomial behaves as 


pfi’ 2 \ cost ) 


Q( A" 1 ) 

t 


uniformly for A 1 < t < tt — A 1 [IS], Theorem 7.32.2], it stems that the ultras- 
pherical polynomial Ca is uniformly bounded as 


\Ch(jp-q)\ < ^ 


for any points p 1 q in § 3 such that 


distsa (p, q) ^ p and distga (p, —q) > p , 


with a constant C p that only depends on the positive constant p. 

Using the formulas of Section [5] it is now easy to show that for any j and any 
fixed positive radius p we have 

„ C 

ll^jllc°(S 3 \(B(P 3 -,p)UB(-P,-,p)) < -J- 


for large A, with a constant that depends on p (and, of course, on v and <5). If we 
now define 

i 


3=1 

and choose p small enough so that the sets B (Pj,p) UB (—Pj,p) are disjoint for all 
j, the same reasoning that we employed in the proof of Proposition 13.21 shows that 


w\ - 

r\ o vp . 1 ( — ^ 

3 

1 J Va J 


C m + 2 (B) 


< cs 


for all 1 ^ i ^ 3 and 1 < j < i, which plays a role completely analogous to that of 
Proposition 13.21 in the generalized context that we are now considering. The rest of 
the argument remains exactly as in Section [3] so the result follows. □ 
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In particular, this yields the existence of Beltrami fields in the sphere having 
prescribed sets of closed vortex lines and tubes (modulo diffeomorphism) around 
any finite number of points P± ,..., Pi. These lines and tubes are contained in balls 
of radius 1/A. This line of reasoning also allow us to prove an analog of Theorem ll.il 
in any lens space L(p,q): 

Theorem 8.2. Let S be a finite union of (pairwise disjoint, but possibly knotted 
and linked) closed curves and tubes contained in a contractible subset of a three- 
dimensional lens space L(p , q). Then for any large enough even integer A there exists 
a Beltrami field u satisfying the equation curl it = Xu and a diffeomorphism $ of 
L(p , q) such that $(5) is a union of vortex lines and vortex tubes ofu. Furthermore, 
this set is structurally stable. 


Proof. The lens space can be written as 


L(p,q)=S 3 /G, 


where G is a finite isometry group isomorphic to Z p . We can assume that G is 
generated by certain isometry g. Let us now fix a point po £ § 3 and set 

Pj ■= 9 3 ' Po 


for 0 < j < p — 1. If if is a patch of normal geodesic coordinates around po, we 
will also set T, (x) := \H(< 7 _J ’ • x). Notice that if p is odd there are not any points 
in the set {Pj}jZ q that are antipodal to each other, while for p even Pj and Pk are 
antipodal if and only if |j — k | = §■ 

Let us fix a Beltrami field v in R 3 as in Theorem 12.21 Theorem 18.11 then ensures 
the existence of a Beltrami field u in § 3 such that 




- Vi 


< e, 


C m (B) 


where 0 ^ j Si p' — 1 with p' := p if p is odd and p' := | if p is even. Here 
vq := v and Vj := 0 for 1 ^ j ^ p' — 1. Notice that, as A is even, we saw in the 
previous remark that u is odd, i.e., u{x) = —u(—x), so that u is equivariant under 
the isometry x i— > —x. Hence, by construction, the vector field 


p'-i 

u '■= 

j=0 


is G- equivariant, and therefore it defines a vector field in the quotient space L(p , q) = 
§ 3 /G that we still denote by u with some abuse of notation. Arguing exactly as 
in the proof of the main theorem one can show that the vector field u on L(p , q) 
indeed has the desired properties, so the statement then follows. □ 


In the torus, the distribution of rational points on the 2-sphere is key. The proof 
that we have given holds provided that the eigenvalue A is an odd integer of suffi¬ 
ciently large absolute value. It does not say anything about even integers, or about 
eigenvalues that are not integers. This assertion can be refined a little, however. 
We have seen that for any eigenvalue A of the curl operator in T 3 there is a set 
of points {£«}(/=! lying on the unit sphere § 2 of R 3 such that A£ n £ 1? (this is 
obvious from the fact that one can write A = |fc| with k £ Z 3 ). Therefore, in the 
proof of Theorem 12.II for the torus (cf. Section [ 7 ]) one can substitute the collection 
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of odd integers A by any subset of eigenvalues A for which there is a set of points 
{£ n }^=i C § 2 (depending on A and such that the rescalings A£ ra are in Z 3 ) that 
becomes dense in the sphere as |A| —> oo along this subset of eigenvalues. In partic¬ 
ular, replacing the density condition by the more stringent assumption that {£„} 
becomes equidistributed on the sphere, it turns out that the characterization of the 
numbers A that satisfy this property is somehow related to the celebrated Linnik 
problem in number theory. In particular, since the aforementioned equidistribution 
property holds for any eigenvalue for which the integer A 2 is square-free [5], we im¬ 
mediately infer that the statement of Theorem 11.11 also holds for any large enough 
eigenvalue A of curl (possibly even or non-integer) for which A 2 is square-free. 
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